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Abstract — We investigate tlie propagation of surface waves 
along a spatially dispersive graphene sheet, including substrate 
effects. The proposed analysis derives the admittances of an 
equivalent circuit of graphene able to handle spatial dispersion, 
using a non-local model of graphene conductivity. Similar to 
frequency selective surfaces, the analytical admittances depend 
on the propagation constant of the waves traveling along the 
sheet. Dispersion relations for the supported TE and TM modes 
are then obtained by applying a transverse resonance equation. 
Application of the method demonstrates that spatial dispersion 
can dramatically affect the propagation of surface plasmons, 
notably modifying their mode confinement and increasing losses, 
even at frequencies where intraband transitions are the dominant 
contribution to graphene conductivity. These results show the 
need for correctly assessing spatial dispersion effects in the 
development of plasmonic devices at the low THz band. 

Index Terms — Surface waves, graphene, spatial dispersion, 
plasmonics 



I. Introduction 

The propagation of electromagnetic waves along graphene 
sheets have recently attracted significant attention IT], 
131 . Graphene, thanks to its interesting electrical and optical 
properties ||4], provides new possibilities for surface wave 
propagation at optical, infrared, and low THz frequencies |l3l, 
0. Specifically, graphene enables the development of novel 
plasmonic devices at relatively low frequencies ||3|, ||6l, Q, 
IS) in contrast with noble metals which only allows plasmonic 
propagation in the visible range ||9l, ifTO) . In addition, graphene 
is inherently tunable, via the application of a electrostatic 
or magnetostatic bias field, leading to novel reconfiguration 
possibilities ||5l, ifTTl . Many research groups have theoretically 
studied the propagation of plasmons (i.e., electromagnetic 
waves propagating at the interface between a conductor and 
a dielecti-ic) along graphene IT], 12, 13, HI], O- Also, 
different configurations have recently been proposed to im- 
prove the characteristics of this propagation, including parallel 
plate pairs lfT2l . IT4l or waveguides [15] . Furthermore, recent 
studies have demonstrated enhanced transmission through a 
stack of monolayer graphene sheets IT6l . However, though 
graphene is known to be spatially dispersive ITT] , spatial 
dispersion has usually been neglected in works related to wave 
propagation at low THz frequencies. In 151, a non-local model 
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of the graphene tensor conductivity, valid in the frequency 
range where intraband contributions of graphene dominate, 
was proposed. Based on this model, l2| also predicted that 
spatial dispersion would become significant for extremely slow 
surface waves. In case of graphene in free space, this occurs 
when interband contributions of graphene dominate (moderate 
to high THz frequencies). 

In this work, we investigate the propagation of surface 
waves along a spatially dispersive graphene sheet, including 
substrate effects, and we show that spatial dispersion affects 
the behavior of these waves even in the low THz range. The 
sheet is characterized using the graphene tensorial conductivity 
obtained by applying the non-local spatially dispersive model 
of graphene presented in 12], valid in the absence of external 
magnetostatic bias field and at frequencies where intraband 
contributions of graphene dominate. The operator components 
of the conductivity tensor are then mapped onto the admit- 
tances of a rigorous Green's function-based equivalent circuit 
of graphene ITSl . These equivalent admittances, which are 
similar to those found in the analysis of frequency selec- 
tive surfaces |fT9] . Il20l . analytically show that the influence 
of spatial dispersion directly depends on the square of the 
wave propagation constant (kp). Then, a transverse resonance 
equation (TRE) 1211 is imposed to compute the dispersion 
relation of surface waves along spatially dispersive graphene. 
Similar to the case of non spatially dispersive graphene sheets, 
transverse electric (TE) and transverse magnetic (TM) modes 
are supported. Analytical dispersion relations are provided for 
TM and TE surface waves propagating along a graphene sheet 
embedded into an homogeneous medium, while approximate 
expressions are given for the case of graphene surrounded by 
two different dielectrics. The derived dispersion relations ana- 
lytically show that the permittivity of the surrounding media, 
operation frequency and spatial dispersion similarly contribute 
to determine the characteristics of kp, suggesting that the 
influence of graphene spatial dispersion in the propagating 
waves may be strongly affected by the environment of the 
sheet. Numerical results confirm that spatial dispersion is an 
important mechanism for wave propagation along graphene 
sheets, leading to surface modes that can significantly differ 
from those found neglecting spatial dispersion. These features, 
which include variations in the mode confinement and higher 
losses, should be rigorously taken into account in the devel- 
opment of novel plasmonic devices at the low THz band. 

The paper is organized as follows. Section|lI]derives the ana- 
lytical relations between the admittances of a Green's function- 
based equivalent circuit of graphene and the components of 
the spatially dispersive conductivity tensor Then, Section HHl 
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Fig. 1: Spatially dispersive graphene sheet on a dielectric 
substrate under arbitrary plane-wave incidence (a) and its 
equivalent transmission line model ifTSl (b). 



computes the dispersion relation of surface waves propagating 
along graphene, providing analytical expressions for the case 
of a graphene sheet embedded into an homogeneous media. 
Section |IV] discusses the characteristics of surface waves along 
spatially dispersive graphene, taking into account the sur- 
rounding media. Finally, conclusions and remarks are provided 
in Section IV] 

II. Equivalent Circuit of a 
Spatially Dispersive Graphene Sheet 

Let us consider an infinitesimally thin graphene sheet in 
the plane z = and separating two media, as illustrated in 
Fig. [Ta] The sheet is characterized by the conductivity tensor 
CT, obtained by applying a spatially dispersive (non-local) 
model of graphene Jll in the absence of external magnetostatic 
biasing fields {Bq = 0). This anisotropic conductivity reads 



Cr(w,/ic,T,r) 



(7x' y' 
<Jyiy' 
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where T is temperature, t is the electron relaxation time, 
yUc is the chemical potential and u) is the angular frequency. 
Due to the spatial dispersion of graphene, the conductivity 
components become operators JS], ll22l 
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fcs is the Boltzmann constant, q,, is the electron charge, wj? is 
the Fermi velocity (« 10^ m/s in graphene), and the subscripts 
"lo" and "sd" have been included to denote local and spatially- 
dispersive (non-local) terms. Note the presence of a minus 
sign in asd, which was not included in ||2| due to a typo 1221 . 
It is worth mentioning that this model is insensitive to the 
orientation of the graphene lattice, and that anisotropy arises 
as the response of the material to the excitation, as usually 
occurs with spatial dispersion JS], ||221 . Consequently, the 
coordinate system employed in the anisotropic conductivity 
{x' — y', see Eqs. ©-(Illi] is related to the excitation, i.e. a^'x' 
and (Jyiyi provides the response parallel to the E^.' and Eyi 
fields, being x' and y' arbitrary perpendicular directions in the 
infinite sheet. Importantly, Eqs. (|2]i-(|6]l are only strictly valid 
at frequencies where the intraband contributions of graphene 
conductivity dominates, usually the low THz regime. Intraband 
contributions correspond to electron transitions between differ- 
ent energy levels in the same band (valence or conduction), 
while interband contributions are related to the absorption of 
an electromagnetic wave by an electron in a state below the 
Fermi level and the subsequent transition to a state in a higher 
band. The latter phenomenon becomes significant as frequency 
increases due to the higher photon energy ||231 . 

A rigorous equivalent circuit of an anisotropic graphene 
sheet sandwiched between two media is shown in Fig. [Tbl 
iflSl . This four-port circuit relates input and output TE and 
TM waves through equivalent shunt admittances {Yj^ and 
YJ'^^ , respectively), and the cross-coupling between the two 
polarizations through voltage-controlled current generators, 
with coefficients yj^/^^^ and yj^^/^^ for TE-TM and TM- 
TE coupling, respectively. The relation between the anisotropic 
conductivity of graphene and these admittances was provided 
in lITSl for the case of a local model of graphene, i.e. when 
anisotropy is due to an external magnetostatic bias field. 
However, spatial dispersion effects have not been considered 
so far 

Here, we analytically obtain the admittances of the equiv- 
alent circuit shown in Fig. [lb] using a non-local (spatially 
dispersive) model of graphene. For this purpose, we include 
the operator components of the tensor conductivity, Eqs. (|2j- 
(|4]i, into the transmission-line network formalism for comput- 
ing dyadic Green's functions in stratified media ifTSl . 11241 . 
For the admittances derivation, we consider a uniform plane 
wave impinging on the graphene sheet from an arbitrary 
direction {9, </>) of medium 1 (see Fig. [Tall, and then we 
use the corresponding auxiliary coordinate system for the 
conductivity tensor. The incoming wave has a wavenumber 
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k = kj-Cx + kyCy + kzCz, where k^ ~ ki cos^cf)) sm{9), 
ky = fci sin((/)) sin(6'), kz = ki cos{6), e^, ey, iz are unit vec- 
tors, and ko = 27r/Ao and ki = y^Eriko are the wavenumbers 
of free space and medium 1, respectively. Imposing boundary 
conditions in the graphene plane and following the approach 
described in ifTSl . 1241 these admittances can be written as 



Yj''ikp)=aio + k-'p[a,d + (3sd] 
(kp) = aio + fc^[a,sd + I3sd] 



TM 



Yr/^'\kp)=Yr'/^^{kp)=0, 



(7) 
(8) 
(9) 



where kp = k'^ + ky is the propagation constant of the wave 
traveling along the graphene sheet. 

The importance of Eqs. (|7]i-(|9]l is two fold. First, they 
analytically indicate that the influence of spatial dispersion 
is directly proportional to k^. Consequently, this phenomenon 
will be more important for very slow waves {kp ko), which 
usually appear at moderately high THz frequencies (where 
interband contributions of graphene dominate and Eqs. (l2)-(|4li 
are not strictly valid |j2l). However, note that slow waves can 
also be obtained in graphene sheets at low THz frequencies by 
using substrates with high permittivity constant (3l- Second, 
the equivalent admittances do not depend on the direction 
of propagation of the wave along the sheet (0) and there 
is no coupling between the TM and TE modes. This is in 
agreement with the non-local model of graphene employed in 
this development ID, 1221 which assumes an isotropic infinite 
surface where spatial dispersion arises as a response to a given 
excitation. 

III. Dispersion Relation for a 
Spatially Dispersive Graphene Sheet 

The dispersion relation of surface waves on a spatially 
dispersive graphene sheet sandwiched between two different 
media can be obtained by imposing a transverse resonance 
equation l2Tl to the equivalent circuit shown in Fig. [lb] The 
solution of the TRE provides the propagation constant kp of 
the surface wave propagating along the sheet. Following the 
approach described in ifTSI . the desired dispersion relation can 
be obtained as 
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are the TE and TM admittances, transverse propagation con- 
stant and wavenumber of media 1 and 2, respectively. 

Eq. ( [Tol l does not generally admit any analytical solution 
and must be solved using numerical methods [25 1. Impor- 
tantly, the mathematical solutions of this equation must be 
carefully checked to verify that they correspond to physical 
modes. Specifically, physical modes must fulfill the law of 
energy conservation, i.e. surface waves cannot be amplified 
while propagating along the structure, and the Sommerfeld 



boundary radiation condition 11211 . Also, note that the solution 
of Eq. (fTOl i leads to the propagation constant of a transverse 
electric (TE) or a transverse magnetic (TM) mode, as in the 
case of isotropic graphene |[T]. These cases are examined 
below. 

A. TM modes 

The dispersion relation for TM modes propagating along a 
spatially dispersive graphene sheet can be obtained by solving 



be expressed as 
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[see Eq. (fT0]i1. This equation may 
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If the graphene sheet is embedded into an homogeneous 
host medium, i.e. e^i = £r2 = £r, the dispersion relation can 
be simplified to 
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which has 6 complex roots but can be solved analytically 
using standard techniques [|26l . Also, considering two different 
media and assuming the usual non-retarded regime {kp ^ ko), 
Eq. (fT2] i reduces to 



kp+ kp- 



-JUJ- 



= 0, 



(14) 



Ctsd + Psd Oisd + Psd 

which is a cubic equation with three complex roots. Note 
that the operation frequency (uj) and the permittivity of the 
surrounding media (e^i and £^2) only appear in the free term 
of Eq. ( fT4] i. where they multiply each other Therefore, the 
behavior of surface waves propagating along a graphene sheet 
is mainly determined by this product, suggesting that similar 
responses will be obtained for larger frequencies if the permit- 
tivity is simultaneously reduced, or vice-versa [Similar conclu- 
sions are reached by closely examining Eq. (fTjil. However, 
this does not mean that frequency and permittivity are fully in- 
terchangeable since graphene conductivity itself is frequency- 
dependent. In the asymptotic limit, ujeo{Sri +£7-2) >> kpaio, 
Eq. (fl4l i can be solved analytically and yields 



(£rl + £r2) 



/ J^^£o- 



asd + Psd 



(15) 



This equation shows that in the limit of very high fre- 
quencies, or a large permittivity of the surrounding media, 
spatial dispersion will be the main phenomenon governing 
wave propagation. However, note that interband contributions 
of graphene conductivity, which are dominant at such high 
frequencies, have not been considered in this analysis. 

In the absence of spatial dispersion (asd = Psd = 0), Eq. 
(fl4l i simpUfies to the usual expression fS] 

(£rl + £r2) 
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Fig. 2: Characteristics of surface waves propagating along a spatially dispersive (SD) graphene sheet versus frequency computed 
using Eq. ( flJl l. Reference results, related to surface waves propagating along non spatially dispersive graphene (|T] computed 
taking into account intraband and intraband+interband contributions of conductivity, are included for comparison purposes, (a) 
and (b) show the normalized propagation constant and losses of surface waves along a free-space standing graphene sheet, 
(c) and (d) show the normalized propagation constant and losses of surface waves along a graphene sheet embedded into an 
homogeneous media with Sr = 11.9. Graphene parameters are /ic = 0.05 eV, r = 0.135 ps and T = 300° K. 



The cut-off frequency of the propagating TM modes can 
be obtained by numerically finding the lowest frequencies 
which fulfill Eq. ( fT2] i. A very good approximation is obtained 
by identifying the frequency range of /?ti((T/o) < 0, which 
is the condition for TM mode propagation along a non- 
spatially dispersive graphene sheet |1)|. Note that the cut-off 
frequency of the supported TM modes mainly depends on the 
characteristics of the local conductivity of graphene, which 
can be externally controlled by tuning the chemical potential 
of graphene via the field effect. 



B. TE modes 

The dispersion relation for TE modes propagating along a 
spatially dispersive graphene sheet can be expressed as 



+ = - [mo + K{asd + Psd)] ■ 

(18) 

In case the graphene sheet is embedded into an homogeneous 
medium, with e^i = £r2 = £r, this equation admits the 
analytical solution shown in Eq. (fT7] i. 

Similarly as in the case of TM modes, the cut off frequency 
of TE modes mainly depends on the local conductivity of 
graphene, aio- Specifically, an accurate condition for the prop- 
agation of these modes along a graphene sheet is Im{aio) > 

m. 
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Fig. 3: Normalized phase constant of surface waves propagating on a spatially dispersive graphene sheet versus frequency, 
computed for different values of (a) chemical potential [using t = 0.135 ps] and (b) relaxation time [using /i^ ~ 0.05 eV]. 
Reference results, related to surface waves on non spatially dispersive graphene sheet computed only taking into account 
intraband contributions of conductivity [[I], are included for comparison purposes. The permittivity of the surrounding media 
is £r = 11.9 and temperature is T = 300° K. 



IV. Numerical Results 

In this section, we investigate the influence of spatial disper- 
sion in the characteristics of surface waves propagating along 
a graphene sheet, taking into account the surrounding media. 
Specifically, we study the normalized propagation constant 
[mode confinement, Re{kp/ko), and losses, Im{kp/kQ)] of 
TM waves along spatially dispersive graphene, and compare 
the results with the ones obtained neglecting spatial dispersion 
lH]. In the numerical study, we focus our results on TM waves, 
which are known to be of interest for plasmonic devices 
(l3|, 161. Though TE surface waves are also supported by 
graphene sheets, they present similar characteristics to waves 
propagating in free space (i.e. kp sa fcg) lH] and thus have 
less practical interest. We compute reference results for non 
spatially dispersive graphene sheet considering intraband and 
intraband+interband contributions of conductivity. The aim 
of this comparison is to identify which phenomenon (spatial 
dispersion or interband contributions) becomes dominant as 
frequency increases. Furthermore, we will also investigate 
the influence of spatial dispersion in the propagating surface 
waves as a function of different parameters of graphene. Our 
numerical simulations consider graphene at T = 300° K and 
a relaxation time of 0.135 ps, in agreement with measured 
values Is), El- The results shown here have been obtained 
numerically [solving Eq. (fT2] il or analytically [from Eq. (fT3]l1. 
depending on the surrounding media of graphene. Besides, 
note that the analytical solution of Eq. (fT4] i leads to results 
with differences smaller than 0.1% with respect to those 
numerically obtained from Eq. (fT2] i. further confirming the 
validity of this expression. 

First, we consider the case of a graphene sheet with chem- 
ical potential = 0.05 eV surrounded by air. Figs. |2a} 
l2b] show the normalized propagation constant and losses of 



surface waves propagating along the spatially dispersive sheet. 
Specifically, two modes are supported by the structure, which 
are the physical solutions of Eq. (fTsT l. The first mode (denoted 
as "SD - mode 1") presents extremely similar characteristics 
as compared to a TM surface mode supported by a non 
spatially dispersive graphene sheet, computed considering only 
intraband contributions of graphene. Figs. |2all2b] also include 
similar computations but considering interband contributions 
as well. The effect of interband contribution become visible 
at high frequencies, increasing the losses of the mode. Thus, 
we can conclude that for this particular case the influence 
of spatial dispersion is negligible in the low THz range. In 
|j2] similar conclusions were obtained analyzing, in the trans- 
formed Fourier domain, the ratio between graphene conduc- 
tivity and the terms related to spatial dispersion. Interestingly, 
here we also observe that a second mode, denoted as "SD - 
mode 2", is supported due to the presence of spatial dispersion. 
It is observed that this mode is extremely lossy in the band 
of interest, which greatly limits its possible use in practical 
applications. 

We present in Figs. |2c]|2dl a similar study, but considering 
now a graphene sheet embedded into an homogeneous material 
with Er ~ 11.9. Results demonstrate that spatial dispersion 
becomes the dominant mechanism of wave propagation in this 
case, even at relatively low frequencies, drastically changing 
the behavior of the two modes supported by the graphene 
sheet. At low frequencies, the first mode ("SD - mode 1") 
presents very similar characteristics as compared to a TM 
surface mode on along a graphene sheet neglecting spatial 
dispersion. However, the confinement and losses of the mode 
increases with frequency. At around 2 THz, losses exponen- 
tially decreases thus leading to a non-physical improper mode 
at frequencies above 2.5 THz. The behavior of the second 
mode ("SD - mode 2") is also affected by spatial dispersion. 
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Fig. 4: Characteristics of surface waves propagating along a spatially dispersive (SD) graphene sheet as a function of surrounding 
media permittivity, computed using Eq. ( fT3] l. Reference results, related to surface waves propagating along non spatially 
dispersive graphene O] computed taking into account intraband and intraband+interband contributions of conductivity, are 
included for comparison purposes, (a) and (b) show the normalized propagation constant and losses of the surface waves at 
1 THz. (c) and (d) show the normalized propagation constant and losses of the surface waves at 3 THz. Graphene parameters 
are /i^ = 0.05 eV, t = 0.135 ps and T = 300° K. 



As frequency increases, the confinement and losses of the 
mode decreases. Above 2 THz, the confinement of the mode 
remains relatively constant and losses slightly increase. It 
should be noted that although the phase constant of modes 
1 and 2 intersect at around 2 THz, their attenuation constants 
are different at that frequency thus allowing to clearly identify 
the modes. Besides, note that the frequency where the phase 
constant of the modes intersect defines the frequency region 
where the influence of spatial dispersion starts to be dominant. 
Furthermore, mode coupling is occurring between the two 
modes supported by the graphene sheet. This coupling is 
governed by the occurrence of a complex-frequency-plane 
branch point that migrates across the real frequency axis ll27l . 

EH). 

Importantly, the influence of spatial dispersion on the char- 
acteristics of the supported modes strongly depends on the 
features of graphene. As an illustration. Fig. |3] presents the 
phase constant of the two modes propagating along a graphene 
sheet (with the parameters employed in Fig. |2c) as a function 
of the chemical potential jic and the relaxation time r. We 
observe in Fig. |3a] that the chemical potential controls the 



frequency where the phase constant of the two modes intersect, 
which is the frequency where spatial dispersion starts to 
have a dominant effect on the behavior of the modes. Note 
that an increase of fic up-shifts this frequency. In addition. 
Fig. [3b] shows that the relaxation time of graphene determines 
the phase constant of the surface waves in the frequency 
region where spatial dispersion is dominant. Increasing t 
modifies the phase constant of the modes, which tend to 
have a similar behavior versus frequency, and reduces their 
attenuation losses. 

For the sake of completeness, we present in Fig. |4]a study of 
surface waves propagating along spatially dispersive graphene, 
but considering now a constant frequency and varying the 
permittivity of the surrounding media. It is worth mentioning 
that in practice the material permittivity might affect graphene 
relaxation time ll23l , but this effect is neglected here for 
convenience. Figs. l4aBbl present the normalized phase and 
attenuation constants of surface waves propagating along 
graphene for a fixed frequency of 1 THz. The first mode 
is extremely similar to a TM surface mode along graphene 
neglecting spatial dispersion. As expected from Eq. (fl4l i. mode 
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Fig. 5: Characteristics of surface waves propagating along a spatially dispersive (SD) graphene sheet deposited on a silicon 
substrate (e^ = 11-9) versus frequency, computed using Eq. (fT2] i. Reference results, related to surface waves propagating 
along non spatially dispersive graphene HI computed taking into account intraband and intraband+interband contributions of 
conductivity, are included for comparison purposes, (a) and (b) show the normalized propagation constant and losses of the 
surface waves, respectively. Graphene parameters are similar to those of Fig. |2] 



confinement and losses increases when the dielectric constant 
is very high. Interesting, there are many similarities in the 
behavior of the different modes in this situation and in the 
previous example, where the surrounding media was constant 
(with a low permittivity value of = 1) and frequency was 
increased. Moreover, we present in Figs. l4cTl4dl the same study 
but at the higher frequency of 3 THz. It is observed that spa- 
tial dispersion becomes the dominant phenomenon for wave 
propagation. In fact, it leads to propagating modes with similar 
behavior as in the previous example (see Figs. I2cll2dt . where 
we used a fixed large value of permittivity (e^ = 11-9) and 
varied frequency. This example confirms that, neglecting the 
frequency dependence of graphene conductivity, permittivity 
and frequency play a similar role in the characteristics of 
surface waves propagating along spatially dispersive graphene. 

Finally, we investigate in Fig. |5] the propagation of sur- 
face waves in a more realistic scenario, which consists of 
a spatially dispersive graphene sheet deposited on a silicon 
substrate (e^ = 11-9), as illustrated in Fig. [Ta] Very similar 
behavior as compared to the previous examples is obtained. 
Furthermore, these results demonstrate that the variation of the 
surrounding media directly controls the frequency range where 
spatial dispersion is noticeable. Analyzing these variations, we 
can conclude that increasing the permittivity of the media 
surrounding graphene down shift the frequencies where the 
spatial dispersion phenomenon dominates wave propagation. 

V. Conclusions 

The propagation of surface waves along a spatially dis- 
persive graphene sheet has been addressed. The analysis, 
that was based on the transverse resonance method extended 
to handle graphene spatial dispersion, allowed to obtain the 
desired dispersion relation for surface waves propagating along 
a graphene sheet. Our results have demonstrated that spatial 
dispersion is an important mechanism which contributes to the 
propagation of surface waves along graphene sheets and that 



its influence cannot be systematically neglected. In fact, the 
presence of spatial dispersion at frequencies where intraband 
contributions of graphene dominate lead to surface waves with 
very different characteristics from those propagating along 
graphene sheets neglecting spatial dispersion. In addition, 
the frequency region where spatial dispersion is noticeable 
depend on the surrounding media and the chemical potential of 
graphene, while the features of the propagating surfaces waves 
are mainly determined by the relaxation time. These results 
demonstrate that the influence of spatial dispersion on surface 
waves propagating on graphene sheets should be rigorously 
taken into account for the development of novel plasmonic 
devices in the low THz range. 

The study presented here has been based on a theoretical 
model of graphene which characterizes it as an infinitesimally 
thin layer with an associated tensor conductivity. This model 
uses the relaxation time approximation, only takes into ac- 
count intraband contributions of graphene and is valid in the 
absence of external magnetostatic biasing fields. Further work 
is needed in this area to analyze the behavior of surface waves 
along spatially dispersive graphene when these assumptions 
are not satisfied. 

Acknowledgement 

This work was supported by the Swiss National Science 
Foundation (SNSF) under grant 133583 and by the EU FP7 
Marie-Curie lEF grant "Marconi", with ref. 300966. The 
authors wish to thank Prof. G. W. Hanson (University of 
Wisconsin-Milwaukee, USA), Dr Garcia- Vigueras and Dr E. 
Sorolla-Rosario (Ecole Polytechnique Federale de Lausanne, 
Switzerland) for fruitful discussions. 

References 

[1] G. W. Hanson, "Dyadic green's functions and guided surface waves for 
a surface conductivity of grapliene," Journal of Applied Physics, vol. 
103, p. 064302, 2008. 



8 



[2] , "Dyadic green's functions for an anisotropic non-local model of 

biased graphene," IEEE Transactions on Antennas and Propagation, 

vol. 56, no. 3, pp. 747-757, March 2009. 
[3] M. Jablan, H. Buljan, and M. Soljacic, "Plasmonics in graphene at 

infrai'ed frequencies," Physical review B, vol. 80, p. 245435, 2009. 
[4] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. 

Dubonos, I. V. Grigorieva, and A. A. Firsov, "Electric field effect in 

atomically thin carbo filts," Science, vol. 306, pp. 666-669, 2004. 
[5] A. Vakil and N. Engheta, "Transformation optics using graphene," 

Science, vol. 332, pp. 1291-1294, 2011. 
[6] W. L. Barnes, A. Dereux, and T. W. Ebbesen, "Surface plasmon 

subwavelength optics," Nature, vol. 424, pp. 824-830, 2003. 
[7] Y. V. Bludov, M. I. Vasilevskiy, and N. M. R. Peres, "Mechanism 

for graphene-based optoelectronic switches by tuning surface plasmon- 

polaritons in monolayer graphene," Europhysics Letters, vol. 92, p. 

68001, 2010. 

[8] M. Tamagnone, J. S. Gmez-Daz, J. R. Mosig, and J. Perruisseau-Carrier, 

"Reconfigurable thz plasmonic antenna concept using a graphene stack," 

Applied Physics Letters, vol. 101, p. 214102, 2012. 
[9] J. M. Pitarke, V. M. Silkin, E. V. Chulkov, and R M. Echenique, "Theory 

of surface plasmons and surface-plasmon polaritons," Rep. Prog. Phys., 

vol. 70, 2007. 

[10] Y. Wang, E. W. Plummer, and K. Kempa, "Foundations of plasmonics," 
Advances in Physics, vol. 60, pp. 799-898, 2011. 

[11] V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, "Magneto- 
optical conductivity in graphene," Journal of Physics: Condensed 
Matter, vol. 19, no. 2, p. 026222, 2007. [Online]. Available: 
http://sta cks.iop.org/0953-8 984/19/i=2/a=026222. 

[12] E. H. Hwang and J. D. Sarma, "Dielectric function, 
screening, and plasmons in two-dimensional graphene," Phys. 
Rev B, vol. 75, p. 205418, May 2007. [Online]. Available: 
http://link.aps.org/doi/10.1103/PhysRevB.75.205418 

[13] J. S. Gomez-Di'az and J. Perruisseau-Carrier, "Propagation of hybrid 
TM-TE plasmons on magnetically-biased graphene sheets," Journal of 
Applied Physics, vol. 112, p. 124906, 2012. 

[14] G. W. Hanson, "Quasi-transverse electromagnetic modes supported by 
a graphene parallel-plate waveguide," Journal of Applied Physics, vol. 
104, p. 084314, 2008. 

[15] J. Christensen, A. Manjavacas, S. Thongrattanasiri, F. H. L. Koppens, 
and F. J. G. de Abajo, "Graphene plasmon waveguiding and hybridiza- 
tion in individual and paired nanoribbons," ACS Nano, vol. 6, p. 431440, 
2011. 

[16] C. S. R. Kaipa, A. B. Yakovlev, G. W. Hanson, Y. R. Padoora, 
F. Medina, and F. Mesa, "Enhanced transmission with a graphene- 
dielectric microstructure at low-terahertz frequencies," Phys. Rev. B, 
vol. 85, p. 245407, 2012. 

[17] L. A. Falkovsky and A. A. Varlamov, "Space-time dispersion of 
graphene conductivity," European Physical Journal B, vol. 56, pp. 281- 
284, 2007. 

[18] G. Lovat, "Equivalent circuit for electromagnetic interaction and trans- 
mission through graphene sheets," IEEE Transactions on Electromag- 
netic Compatibility, vol. 54, pp. 101-109, February 2012. 

[19] B. A. Munk, Frequency Selective Surfaces: Theory and Design. John 
Wiley and Sons, 2000. 

[20] S. Maci, M. Caiazzo, A. Cucini, and M. Casaletti, "A pole-zero matching 
method for ebg surfaces composed of a dipole fss printed on a grounded 
dielectric slab," IEEE Transactions on Antennas and Propagation, 
vol. 53, pp. 70-81, Januaiy 2005. 

[21] R. E. Collin and F. J. Zucker, Antenna Theory. McGraw-Hill, 1969. 

[22] G. W. Hanson, Private Communication, 2012. 

[23] M. Dressel and G. Gruner, Electrodynamics of Solids. Cambridge, 
U.K.: Cambridge Univ. Press, 2002. 

[24] K. A. Michalski and J. R. Mosig, "Multilayered media Green's functions 
in integral equation formulations," IEEE Transactions on Antennas and 
Propagation, vol. 45, no. 3, pp. 508-519, March 1997. 

[25] W. H. Press, S. A. Teukolsky, W. T. Vetteriing, and B. R Flannery, Nu- 
merical Recipes in Fortran 90, The Art of Parallel Scientific Computing. 
Cambridge University Press, 1996. 

[26] L. A. Pipes and L. R. Harvill, Applied Mathematics for Engineers and 
Physicist, 3rd ed. McGraw Hill, 1971. 

[27] G. W. Hanson and A. B. Yakovlev, "Investigation of mode interaction on 
planar dielectric waveguides with loss and gain," Radio Science, vol. 34, 
pp. 1349-1359, 1999." 

[28] A. B. Yakovlev and G. W. Hanson, "Mode transformation and mode 
continuation regimes on waveguiding structures," IEEE Transactions on 
Microwave Theory and Techniques, vol. 48, pp. 61-15, 2000. 



